Antibaryon density in the central rapidity region of a heavy-ion collision 
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We consider (anti-)baryons production in heavy ion col- 
lisions as production of topological defects during the chi- 
ral phase transition. Non-zero quark masses which explic- 
itly break chiral symmetry supress the (anti-)baryon density. 
Hardly any (anti-)baryons will be produced in the central ra- 
pidity region of a heavy ion collision. 

In the series of papers Ellis and Kowalski |JjJ and El- 
lis, Kowalski and Heinz Q considered baryon and an- 
tibaryon production using a topological model inspired 
by the Skyrme picture of baryons as topological defects 
in the quark-antiquark condensate HU . The key idea is 
that the hadronization process is associated with sponta- 
neous chiral symmetry breaking. During these processes 
domains are formed in which the order parameter, the 
quark condensate (0\qq\0), takes the vacuum expectation 
value in each domain independently. Defects may form if 
orientation of the order parameter in adjacent domains 
has a topologically non-trivial configuration. These de- 
fects are interpreted as baryons. Geometrical considera- 
tion in the case of heavy-ion collisions predicts 0.17-0.35 
baryons/fm 3 ||. 

In this communication we would like to consider the 
same scenario under the assumption that the chiral phase 
transition is described by the 0(4) sigma model nearby 
the critical temperature f^d^f. The massless two flavours 
case was considered by Gill |]8| who used Zurek scenario 
H to estimate the density of topological defects after a 
quench. However, it is known fl(i[| that explicit chiral 
symmetry breaking through the non-zero quark masses 
leads to an important modification of the Zurek picture. 
We will show that this effect results in substantial sup- 
pression of the baryon production. 

In the Bjorken scenario JTTJ ] the quark gluon plasma 
comes to thermal equilibrium at the temperature Tj 
around 77 — lfm after an ultrarelativistic heavy ion col- 
lision. As the ultrarelativistic plasma is expanding, it 
cools down adiabatically according to 



T = T / (-) 1/3 

T 



(1) 



w2 , 



is the proper time. Quarks and gluons are 



the relevant degrees of freedom as long as the tempera- 
ture does not drop below T con f, where the confinement - 
deconfinement phase transition takes place. Even above 
T C onf the chiral fields can be defined and measured with 
appropriate < qq > expectation values. One can mea- 
sure e.g. some pion fluctuations. However, the pions do 
not have any dynamics of their own, they are mere phan- 



toms which come out from appropriate quark-antiquark 
correlations in the quark-gluon plasma. As a result above 
T C onf the chiral fields are by definition in thermal equi- 
librium with the thermalized plasma. 

The situation changes at T — T con f. The plasma dis- 
appears and the chiral degrees of freedom become rel- 
evant. If confinement - deconfincmcnt phase transition 
is first order the takeover proceeds at a constant tem- 
perature: latent heat is released and the system cools 
down by adiabatic expansion. The chiral fields enter the 
T < T con f stage in the state of thermal equilibrium and 
they begin to evolve without any quark-gluon heat bath. 

The free evolution of the chiral fields is described by the 
linear 0(4) sigma model with an explicit chiral symmetry 
breaking term 



L = M a d^ a 



A, 



- V 



Ha 



(2) 



where the vector <f> — {a, 7?) is an 0(4) multiplet of real 
scalar fields. The vector 7? represents the pion fields. The 
typical zero temperature fits of some of the parameters 



A = 20 , 

H = (119 MeV) 3 , 
A = A(87.4 MeV) 2 



(3) 



The parameter A depends on temperature. It vanishes 
at T ch = 160-200MeV. If H = 0, T ch would be the tem- 
perature of the second order chiral symmetry breaking 
phase transition. For H ^ there is no phase transition 
but a crossover; <7o(T) =< a >t is nonzero and it varies 
smoothly with temperature, see Fig.l. The pion and 
sigma masses can be extracted from small fluctuations 
around this minimum, see Fig. 2. The masses remain 
substantially nonzero nearby T c h- Correlation lengths 
remain finite and there is no critical slowing down. 

We assume that T con f « T c h- The initial conditions at 
T c h are given by thermal equilibrium distributions. These 
initial conditions themselves rule out any possibility of 
baryon production. 

The effective potential at T c h = 160MeV where 
A(T c h) = is shown in Fig. 3. It has a unique mini- 
mum at ao(T c h) = 44MeV, the sigma mass is M a (T c h) = 
340MeV. The correlation length of the a field is given 
by £cr = l/M a (T c h) « 0.6fm, this is roughly the size 
of the correlated domain of this field. The mass of the 
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pion field is M n (T c h) = 196MeV, its correlation length is 
£ w = l/M w (T c h) = lfm. We define the field a, which is 
an average of a over its correlation domain. In our calcu- 
lations the average is defined by a cut-off in momentum 
space at M a . Fluctuation squared of a around its average 
is given by 



<[*- a Q (T ch )} 2 >-- 



ch 



(2tt) 3 M 2 + k 2 



T ch M a {T ch ) _ TT 
2tt 2 1 4 j 



(4) 



A gaussian distribution f{a;fx) of these domain fluctu- 
ations is shown on an insert in Fig. 3. The same figure 
shows a cross section through the potential at zero tem- 
perature. As the temperature drops below T con f sa T c y L 
there is no quark-gluon plasma to provide an external 
heat-bath for the chiral fields. They evolve freely in 
the zero temperature potential from the initial conditions 
given by thermal equilibrium distributions at T c h- 

From a topological point of view the skyrmion is a map 
from the space compactified to S 3 to the S 3 bottom of 
the zero temperature sombrero potential. As the vector 
field wraps around the target sphere, there must be such 
a point that the field points in the (a, tt) — (—1,0) direc- 
tion. The position of a skyrmion can be identified with 
this point. At T c h there are £ CT -sized domains of a and 
^-sized domains of tt. The field at a given point will 
roll down to the bottom of the sombrero potential in the 
(—1,0) direction, if the domain averaged tt — and at 
the same time the domain averaged a is to the left of the 
top of the zero temperature potential, a < — 11.2MeV. 
The second condition is satisfied with the probability 



P[a < -11.2MeV] 



f (a; fx)da = 0.01, (5) 



where /((t;/i) is the gaussian distribution with fluctua- 
tion given by (|]). The density of the points such that 
the domain averaged tt = can be worked out with the 
general formula |l2l,|l3| 



N\tt 



= 1 r <-i?(x) V 2 tt(x) > -, 3/2 

TT 2 < 7?(x) ff(x) > 



} 



3/2 



M 3 3tt - 8 
TT 2 * 12 - 3tt 



13/2 



0.04 



£ 3 



(6) 



where M n — M^{T c h) and £ ff is the pion correlation 
length at T c h- The ratio of the second derivative of the 
correlation function to the correlation function itself is 
a proper measure of the density of zeros. Indeed if the 
order parameter takes non-zero value and oscillates with 
the small amplitude the correlation function is large in 



compare to its second derivative and as a result its ratio 
is small. In the case when the order parameter oscillates 
near zero with small frequency the correlation function is 
comparable with its second derivative and density of ze- 
ros is of order unity while for large frequency the second 
derivative of the correlation function is larger then the 
function itself and density of zeros is substantial. The 
power 3/2 is |^ for dimensional reasons. 

As our definition does not distinguish baryons from 
antibaryons, the expected density of antibaryons is 



N s = - x P[a < -11.2MeV] x N[tt = 0] = 



= 2-10 4 antibaryons/fm 3 
for the assumed T ch = 160MeV. 



(7) 



Discussion 

We found negligible (anti-) baryons density (of order 
10 _4 antibaryons/fm 3 ) after the chiral symmetry phase 
transition in quark-gluon plasma at zero baryon chemi- 
cal potential. The (anti-)baryon density suppression re- 
sults from the initial condition of the order parameter 
which is in thermal equilibrium at the critical point and 
from the non zero light quark masses which substantially 
bias these initial conditions. The probability of the topo- 
logical defect formation is exponentialy supressed by the 
explicit symmetry breaking parameter present in the free 
energy describing the phase transition. This is the main 
source of discrepancy between our prediction and that of 
Ellis, Kowalski and Heinz ||. 

We assume that the chiral fields are relevant degrees of 
freedom in the vicinity of critical temperature and that 
the transition is a smooth crossover. We expect that 
even in the case of the first order phase transition (result- 
ing from the strange quark mass and/or nonzero baryon 
chemical potential) symmetry breaking bias would be 
strong enough to suppress baryon production. This issue 
deserves a more careful investigation. 

Our result may overestimate the (anti-)baryon density. 
As the chiral fields cool down the skyrmions can unwind 
through the top of the sombrero potential. This process 
lasts until the system cools down to its Ginzburg tem- 
perature and can vastly decimate the defects. This fur- 
ther reduction takes place if unwinding is not forbiden. 
However, according to Holzwarth the unwinding must be 
forbiden for the sake of baryon number conservation. In 
his model |Q the (a, tt) — (0, 0) point is removed from 
the configuration space of the linear sigma model. In this 
case our estimate is accurate. 

The result of our calculation leads us to the conclusion 
that hardly any baryons will be produced in the central 
rapidity region of a heavy ion collision. 
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Fig. 2. Temperature [MeV] dependence of the 
a (top) and 7r (bottom) masses [MeV] under the 
same model assumptions as in Fig.l. 
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Fig.3. The effective potential V e ff(a,0) [MeV 4 ] 
as a function of a [MeV] at T = T c h (top) and at 
T = (bottom) - bold lines. The gaussian distri- 
bution of (7 at T c h - thin line. 
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Fig.l. ero(T) =< a- >t [MeV] as a function of 
temperature [MeV] for the linear a model with 
the effective A eff (T) = A(l - T/T ch ) and T ch 
Kid MeV. 
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